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Abstract. We present here the first-principles dynamical CPA (coherent potential 
approximation) combined with the tight-binding LMTO LDA+U method towards quantitative 
calculations of the electronic structure and magnetism at finite temperatures in transition metals 
and compounds. The theory takes into account the single-site dynamical charge and spin 
fluctuations using the functional integral technique as well as an effective medium. Numerical 
results for Fe, Co, and Ni show that the theory explains quantitatively the high-temperature 
properties such as the effective Bohr magneton numbers and the excitation spectra in the 
paramagnetic state, and describes the Curie temperatures semiquantitatively. 



Quantitative description of the magnetic properties of materials has been one of the goals in 
magnetism. For the ground-state properties of transition metals and alloys, such a calculation 
scheme has been realized by the density functional theory (DFT). The DFT quantitatively 
explains the ground-state magnetizations of Fe, Co, and Ni as well as their cohesive properties. 
At finite temperatures, on the other hand, it does not explain the local-moment behavior such 
as the Curie- Weiss susceptibilities, and overestimates the Curie temperature by a factor of ten 
due to neglect of the spin fluctuations. Because of these difficulties, quantitative description of 
the finite temperature magnetism has been a long-standing problem in the metallic magnetism 

Quite recently, we have proposed the dynamical coherent potential approximation (CPA) 
[2] combined with the first-principles tight-binding (TB) linear muffin-tin orbital (LMTO) 
method [3] with LDA(Local Density Approximation) +U potential [1] for band calculations, 
and presented the numerical results for Fe and Ni taking into account the dynamical corrections 
up to the second order [5]. In this presentation, we extend the calculations up to the 4th order, 
and demonstrate that the theory explains the finite-temperature properties of transition metals 
quantitatively or semiquantitatively. 

The dynamical CPA is a dynamical version of the single-site spin fluctuation theory developed 
by Cyrot, Hubbard, and Hasegawa since early in the 1970, and has recently been shown to be 
equivalent to the dynamical mean field theory (DMFT) [HE]. The theory describes the magnetic 
properties at finite temperatures efficiently taking into account the dynamical corrections to 
the static approximation which is exact in the high temperature limit. Moreover the present 
approach can treat the transverse spin fluctuations for arbitrary d electron number, which is 
different from the early quantum Monte-Carlo (QMC) calculations combined with the DMFT 
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We adopted the Hamiltonian H = H$ + H\. The TB-LMTO Hamiltonian Hq is written as 



iLa iLjL'cr 

Here e® L is an atomic level on site i and orbital L, tiLju is a transfer integral between orbitals 
iL and jL'. L = (l,m) takes s, p, and d orbitals. aL (aiLa) is the creation (annihilation) 
operator for an electron with orbital L and spin a on site i, and nn, a = a\ LtT aiL a is a charge 
density operator for electrons with orbital L and spin a on site i. It should be noted that 
e® L is not identical with the atomic level en in the DFT-LDA potential because the latter 
contains exchange correlation contributions. For the atomic d level, we adopted the formula 
e° iL = en — dE^ DA /dniLo- @]- Here riiLa is the charge density at the ground state, .E^da * s a 
LDA functional to the intraatomic Coulomb interactions. We adopted the Hartree-Fock type 
form for -Elda proposed by Anisimov et al. [8] because we treat here the itinerant electron 
system where the ratio of the Coulomb interaction to the d band width is not larger than one. 
For the other orbital electrons, we did not make any corrections for simplicity. 
The intraatomic Coulomb interactions H\ is given as 

H% = [ Up numinumi + (U\ - -J)hu m hu m > - ^ JSilm-Silm' ■ (2) 

i '<n m>m' ~ m>m' 

Here Uq (U±) and J in the interaction Hi are the intra-orbital (inter-orbital) Coulomb interaction 
and the exchange interaction, respectively. hu m (su m ) with I = 2 is the charge (spin) density 
operator for d electrons on site i and orbital m. 

We have presented the first-principles dynamical CPA in our previous paper [5]. We briefly 
explain here the outline. In the theory, we transform the interacting Hamiltonian Hi into a 
dynamical potential v in the free energy adopting the functional integral method. Introducing a 
uniform medium, (i.e. a coherent potential) ££ CT (iu; n ), we expand the free energy with respect 
to sites. Note that oj n denotes the Matsubara frequency. The first term in the expansion is the 
free energy for a uniform medium, The second term is an impurity contribution to the 

free energy. The dynamical CPA neglects the higher-order terms, so that the free energy per 
atom is given by 

^cpa = ^[E] - /T C j d£ e-^flKO . (3) 

Here (3 is the inverse temperature. C is a normalization constant. £ is the static field variable on 
a site. i£ e ff(£) is an effective potential projected onto the static field £. It consists of the static 
contribution and the dynamical correction term E^ yn (^). We obtain the latter using the 

harmonic approximation (see Ref. [S] for details of expressions of -E'st(^) and -E'dyn(O)- 

The effective medium, i.e., the coherent potential can be determined by the stationary 
condition 8Tcp^./ST, = 0. This yields the CPA equation as 

(GUtiui)} =F La (iu l ) . (4) 

Note that ( ) at the r.h.s. is a classical average taken with respect to the effective 
potential -E e ff(£)> Fi a (iuJi) = [(iui — Hq — E)~ ]iLaiLo- is the coherent Green function, 
where (H o)iLajL' a is the one-electron TB-LMTO Hamiltonian matrix, and {^)iLajL'a = 
^^(iuj^SijS^Li . Furthermore, Gl ct (£, iu>i) is the impurity Green function given by 



(5) 



Figure 1. Magnetization (M), inverse 
susceptibility (x -1 ), and amplitude of local 
moment ((m 2 ) 1 / 2 ) vs. temperature (T) 
curves for Fe. The dynamical results are 
shown by the solid curves, while the results 
in the static approximation are shown by 
dotted curves. The DMFT results without 
transverse spin fluctuations |7J are shown 
by open squares. Experimental data of 
magnetization [10] are shown by +. 

Here the first term at the r.h.s. is the Green function in the static approximation. The second 
term is the dynamical corrections, and KL a (iuJi) = 1 — Fz ja (iuji)~ 2 5Fi j0 -(iijJi)/ST 1 i jrT (iuji). 

Solving the CPA equation (4) self-consistently, we obtain the effective medium. The magnetic 
moment is then given by {m z L ) = J2na a ^L<j{i^n)i and the density of states (DOS) is obtained 
by means of a numerical analytic continuation. 

In the numerical calculations, the average Coulomb and exchange energy param- 
eters (U and J) are taken from the LDA+U band calculations [U [9]; (U,J) = 
(0.169, 0.066), (0.245, 0.069), (0.221, 0.066) Ry for Fe, Co,_and Ni, respectively. U and U x are 
obtained from the relations Uq = U + 8J/5 and U\ = U — 2J/5. These sets of parameters 
yield the Hartree-Fock Curie temperatures: 12200K, 12100K, and 4950K for Fe, Co, and Ni, 
respectively. 

We present in Fig. 1 the results of calculations for Fe. The calculated susceptibility follows 
the Curie-Weiss law with the effective Bohr magneton number m e g = 3.0/Ub in agreement 
with the experimental value 3.2^b- We find the Curie temperature Tq = 2070K in the static 
approximation. By taking into account the dynamical corrections, the Curie temperature reduces 
to 1930 K, which is comparable to 1900K obtained by the QMC calculations of the DMFT 
without transverse exchange interaction. The theoretical Tq is higher than the experimental 
result 1043K by a factor of 1.8. The ground-state magnetization obtained by an extrapolation, 
2.58 fiB is also considerably overestimated as compared with the experimental value 2.22 /^b- 

In the case of Ni, we obtained the Curie temperature Tq = 1420K in the static approximation 
as shown in Fig. 2. The dynamical corrections much reduce the Curie temperature; Tq = 620K, 
which agrees well with the experimental value 630K. Calculated inverse susceptibility shows 
slightly upward convexity. We obtained the effective Bohr magneton number m e g = 1.6/Ub in 
the high temperature region. The result agrees with the experimental one 1.6/xb- As seen from 
Fig. 3, a large reduction of Tq in the dynamical CPA calculation seems to be related to the 
reduction of the DOS p(0) at the Fermi level when we take into account the third and fourth 
order dynamical corrections. The behavior is consistent with the argument of the stability of 
ferromagnetism at the ground state |14j . 

We have also investigated the finite temperature magnetism for the fee Co. In this case, we 
obtained Tq = 2550K which is larger than the experimental value 1388K by a factor of 1.8. 
Calculated effective Bohr magneton number of Co is 3.0/j.b, being in good agreement with the 
experimental value 3.15/Ub- We have also calculated the densities of states from Sc to Cu in the 
high-temperature region, where the present theory works best, and found that the DOS explain 
well the XPS and BIS data. The correlation effects on these excitation spectra will be discussed 
in a separate paper. 

In summary, we have presented the dynamical CPA combined with the first-principles 
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Figure 2. Magnetization, inverse suscepti- 
bility, and amplitude of local moment as a 
function of temperature for Ni. The dynami- 
cal results are shown by the solid curves, while 
the results in the static calculation are shown 
by dotted curves. Dashed curves show the re- 
sults of the 2nd-order dynamical CPA. The 
DMFT results [7] are shown by open squares, 
and experimental data of magnetization [llj 
are shown by +. 
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Figure 3. Calculated DOS in the param- 
agnetic Ni. Solid curve: dynamical CPA, 
dashed curve: 2nd-order dynamical CPA, 
dotted curves: XPS and BIS experimental 
data [121 E3]. 



TB-LMTO LDA+U Hamiltonian. Calculated Curie temperature Tq for Ni agrees with the 
experiment, but Tq for Fe and Co are higher than the experimental ones by a factor of 1.8. On 
the other hand, the present theory explains quantitatively the effective Bohr magneton numbers 
as well as the excitation spectra in the high temperature region. In order to obtain observed Tq, 
we have to take into account the nonlocal correlations definitely. Such an attempt is in progress. 
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Institute of Solid State Physics, University of Tokyo. 
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